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Abstract 

The Maxwell-Dirac equations are the equations for electronic matter, the 
“classical” theory underlying QED. The system combines the Dirac equa¬ 
tions with the Maxwell equations sourced by the Dirac current. 

A stationary Maxwell-Dirac system has p = with <j) independent 

of t. The system is said to be isolated if the dependent variables obey quite 
weak regularity and decay conditions. In this paper we prove the following 
results for isolated, stationary Maxwell-Dirac systems, 

• there are no embedded eigenvalues in the essential spectrum, i.e. —m < 
E <m\ 

• if \E\ < m then the Dirac field decays exponentially as |a:| ^ oo; 

• if |i5| = m then the system is “asymptotically” static and decays ex¬ 
ponentially if the total charge is non-zero. 


1 Introduction 

The Maxwell-Dirac system consists of the Dirac equation 

7“(da — f + = 0, (1) 

with electromagnetic interaction given by the potential Aa', and the Maxwell equa¬ 
tions (sourced by the Dirac current, j“), 

— ^ctAp O^Aq^j 

d°Fai3 = -47rej/3 = -dTreby/jV’- (2) 
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Most studies of the Dirac equation treat the electromagnetic field as given 
and ignore the Dirac current as a source for the Maxwell equations, i.e. these 
treatments ignore the electron “self-field”. A comprehensive survey of these results 
can be found in the book by Thaller Q. This is not surprising, inclusion of the 
electron self-field via the Dirac current leads to a very difficult, highly non-linear 
set of partial differential equations. So difficult in fact that the existence theory and 
solution of the Cauchy problem was not completed until 1997 - seventy years after 
Dirac first wrote down his equation! In a stunning piece of non-linear analysis , 
worked out over an almost twenty year period, Flato, Simon and Tallin (Q) solved 
the Cauchy problem for small initial data. Other contributors to this work on the 
existence of solutions would include Gross [^, Chadam j|], Georgiev §, Esteban 
et iP, and Bournaveas @]. 

There are no known non-trivial, exact solutions to the Maxwell-Dirac equa¬ 
tions in 1 -|- 3 dimensions - all known solutions involve some numerical work. 
These solutions do, however, exhibit interesting non-linear behaviour which would 
not have been apparent through perturbation expansions. The particular solutions 
found in and exhibit just this sort of behaviour - localisation and charge 
screening. See also Das 0 and the more recent work of Finster, Smoller and Yau 

0 - „ 

Finster, Smoller and Yau also point out in that solving the system 

(Einstein-Maxwell-Dirac system in their case) gives, in effect, all the Feynman 
diagrams of the quantum field theory, with the exception of the fermionic loop di¬ 
agrams. Study of the Maxwell-Dirac system should provide an interesting insight 
into non-perturbative QED. 

The aim of the present work is to obtain qualitative information on stationary 
solutions of the Maxwell-Dirac system, in doing so one would hope to be able to say 
something about Maxwell-Dirac models of some simple natural systems involving 
the electron. In fact, it is surprising that so little is understood about these systems. 
There is no Maxwell-Dirac model for a single isolated electron. If we compare the 
situation to that in the other great physical theory of the 20th century, General 
Relativity, the situation could not be more stark. There are a host of solutions to 
Einstein’s equations representing single, isolated gravitating bodies. 

In this article we will prove the following result. 

Main Theorem A stationary, isolated Maxwell-Dirac system has no embedded 
eigenvalues, i.e. —m < E <m. 

If \E\ < m then the Dirac field, decays exponentially as \x\ oo. 

If \E\ = m then the system is “asymptotically static” and ,if the system has non¬ 
vanishing total charge, decays exponentially as |a;| —> oo. 

The content of the theorem will be established through theorems 1 to 5 in 
sections 4 to 7. 

The paper is organised as follows. First, we give a brief overview of 2-spinor 
methods applied to the Maxwell-Dirac equations (some details of 2-spinor calculus 
may also be found in the appendices). In §3 the definitions of stationary, isolated 
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systems are given and some simple consequences are explored. In §4 we address the 
important question of embedded eigenvalues in the spectrum of the Dirac operator, 
the main results being proposition 1 and theorem 1. In §5 we examine regularity 
and decay issues with the main results given in theorems 2 and 3. The next section, 
§6, looks at the special case \E\ = m and its relation to a generalisation of the 
static systems of |^. In §7 we prove the exponential decay of the Dirac field in 
the \E\ = m case. Finally, in §8, we conclude with a brief discussion of the results. 


2 The Maxwell-Dirac Equations 


In this section we give a very brief account of the 2-spinor formulation of the 
Maxwell-Dirac equations, details may be found in Q, some of the “mechanics” of 
the 2-spinor formalism are collected in Appendix A at the end of the paper. 

In [D the 2-spinor form of the Dirac equations was employed to solve (]^) for 
the electromagnetic potential, under the non-degeneracy condition j°‘ja 7 ^ 0. In 
terms of 2-spinors - see below - the non-degeneracy condition can be written as 
ucv^ 7 ^ 0, since j^^ja = 2\ucv^\‘^. Requiring A“ to be a real four-vector gives a 
set of partial differential equations in the Dirac field alone, the reality conditions. 

For 2-spinors ua and (see for an exposition of the 2-spinor formalism) 
we have 


j > with 

ucv'^ 7 ^ 0 (non-degeneracy), 

where A, B = 0,1, A, B — 0,1 are two-spinor indices. The Dirac equations are 

- i e A^^)ua + = 0, 

v2 

(^Ad ^ g A^‘^)va + = 0, 

V2 

where = a°‘^^da, and A^^ = a°‘^"^Aa; here are the Infeld-van der 

Waerden symbols. 

The electromagnetic potential is (see for details), 


(3) 




+ v^v^) 

V2 


e{u‘^Vc) 
The reality conditions are, 


(4) 


d^^iuAUA) = - «%) (5) 

d^^{vAVA) = ^{u^vc -U^Vfs) 

UAd^^VA — VA^^^UA = 0 . 
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The Maxwell equations are, 


The equations (|5) and 
Dirac equations, M and (|). 


= -47rej> = -ATTeV2ap^{uAUx + vav^). (6) 

are entirely equivalent to the original Maxwell- 


3 Isolated, Stationary Maxwell-Dirac Systems 

We recall the definitions of |Q for stationary and isolated systems. 

Definition 1 A Maxwell-Dirac system is said to be stationary if there is a gauge 
in which if = with the bi-spinor f independent of t. Such a gauge will be 

referred to as a stationary gauge. 

Clearly, a stationary gauge is not unique - any gauge transformation if 
6*“*'!/' leaves the system in a stationary gauge. Note that under such a gauge change 
we have, A°‘ —> A°‘ — jSq. We are interested in isolated systems, i.e. systems for 
which the fields decay suitably as |a:| ^ oo, in this case we will require that 
^ 0 as |a;| ^ oo in some stationary gauge. In this particular gauge we will 
write if = for the stationary gauge in which y4° ^ 0 as \x\ oo. Notice 

that for any stationary system in a stationary gauge A°‘ is independent of time, t, 
(see equation (^)). 

In most physical processes that we would wish to model using the Maxwell- 
Dirac system we would be interested in isolated systems - systems where the fields 
and sources are largely confined to a compact region of . This requires that the 
fields decay sufficiently quickly as \x\ oo. 

The best language for the discussion of such decay conditions and other 
regularity issues is the language of weighted function spaces; specifically weighted 
classical and Sobolev spaces. In Q the weighted Sobolev spaces, Wg’^, were used 
following the definitions of [^. These definitions have the advantage that the 
decay rate is explicit: under appropriate circumstances a function in behaves 
as |a;|^ for large |a;|. An element, /, of Wg'^ has in LP for each 

multi-index a for which 0 < |q;| < k] here cr = g/l -\- |a;p and we are working on 
(or some appropriate subset thereof) - see Q or Q and Q (the later papers 
use a different indexing of the Sobolev spaces). We will make use of the Sobolev 
inequality and frequent use of the multiplication lemma. 

Sobolev Inequality (see 1^, [0) If f e Wg'^ then 

(i) ||/||^_^,5 < C'll/llfc,,,^, if n-kp>0 and p < q < 

(ii) ||/||oo.5 < C'll/IU.p.5, ifn-kp<0, and 
|/(a:)| = o(r‘^) as r ^ oo. 
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Multiplication Lemma (see |p^ ) Pointwise multiplication on Ep is a continuous 
bilinear mapping 

if ki,k 2 > k, k < ki + k 2 — nj^, and (5 > + (52- 

We will be interested in the asymptotic region (spatially) of the Maxwell- 
Dirac system, which we denote by Ep = M.^\Bp, where Bp is the ball of radius p. 
We will take our fields to be elements of the function spaces Wg’^{Ep) for certain 
values of the indices k, p and 6. Before introducing the precise definition of an 
isolated system we must (following |l4|) introduce some notation. 

Suppose we have a stationary system and we are in a stationary gauge for 
which ^ 0 as |a:| —> 00 . Write, ua = e~^^*UA and with Ua, 

Va and A°‘ all independent of time, t. Note that ucv^ = UcV^ is a gauge and 
Lorentz invariant complex scalar function, this means we can introduce a (unique 
up to sign) “spinor dyad” {o^, is} with i^oa = 1 - some facts on 2-spinor dyads 
are collected in Appendix A at the end of the paper. The dyad is defined as follows, 
let UcV^ = i?e*x - where R and x are real functions - then write, 

Ua = VRe^^OA and Va = \/RU^ca- 


Note that we must have R > 0 (almost everywhere) because of our non-degeneracy 
condition. 

We can now define our isolated systems, note this definition is a little more 
general than the definition of |Q. 

Definition 2 A stationary Maxwell-Dirac system will be said to be isolated if, in 
some stationary gauge, we have 


if = 


— lEt 


Vr 


e 2 OA 

— i2L-A 

e L 


with E constant and \/R € Wf,’^(Ep); e^^OA, c^^la & 
A“ G Wf'ij^^{Ep), for some t > p > 0 and any e > 0. 

Remarks 


• This definition ensures, after use of the Sobolev inequality and the multipli¬ 
cation lemma, that ip = o(r“'”+'^) and A“ = o(r“^+'^). 

• Notice our condition places regularity restrictions on the fields in the region 
Ep only. In the interior of Bp there are no regularity assumptions. 

• A minimal condition that one may impose on the Dirac field is that it have 
finite total charge in the region Ep, this amounts to 

/ fdx = [ {\Uo\^ + \Ui\^ + \V°\^ + \V^\^) dx 

JEp JEp 

= [ R{\oofdx < 00 . 

JEp 
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This, of course, simply means that Ua and are in L'^{Ep). So Ua and 
would have decay at infinity; roughly, they would decay faster than 
|a;|“ 2 , i.e. we require at least r > |. 

• The spherically symmetric solution of ||] provides an excellent example of 
an isolated, stationary and static Maxwell-Dirac system. 


With our assumption that the Maxwell-Dirac system is isolated and sta¬ 
tionary we can impose the Lorenz gauge condition, without altering the above 
regularity and decay assumptions on To see this we note that 






i=i 


and that the Laplacian, A, gives an isomorphism W^ 2 +e >^^7 assume 

e < 1) - see |^, (and also 1^, |]^ and @). This means there is a unique 
solution, n G W^^‘^{Ep), of the equation 


An -K = 0. 

Consequently, for the gauge change Aa ^ Aa = + da^ we still have G 

W-i+e , with satisfying the Lorenz gauge condition. 

The electromagnetic potential of our stationary, isolated Maxwell-Dirac sys¬ 
tem will be taken to satisfy the two equations. 


Ayl“ 




To end this section 
following sections. 


= 4716^20-“^^i? (OAO^ + i^ai^a) 

A dA^ ^ 




(7) 

( 8 ) 


we present a simple result which we will need in the 


Lemma 1 For a stationary and isolated Maxwell-Dirac system, in the Lorenz 
gauge, 


A^ — G W^p(Ep), go a constant and 


-T] \^p 


A^ G W^’,;(Ep), j = 1,2,3; and t? = 2(r - 1) > 1. 


(9) 

( 10 ) 


Proof. Firstly we note that the source term of the Maxwell equation ( 0 ) is in 
VVi^ 2 .,._^. 2 e{Ep), i.e. G Wi^ 2 .,._^. 2 e{Ep). The Laplacian gives an isomorphism be¬ 
tween and W^ 2 T-\- 2 e see Q, @ (also [^, and @). So there 

exists an a“ G such that, 

Aa“ = 47rej“. 
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Now, as G we have 

+ a“, 

|x| 

where the qa are constants. Applying the Lorenz condition we find qj = 0, for 
j = 1,2,3. □ 

We can improve the decay rates here to higher (negative) order harmonic 
polynomials, at the expense of regularity, by using the Dirac equations to get 
G - but lemma 1^ is sufficient for our purposes. 

The constant qo is the total electric charge of the system (i.e. the electric 
charge of the Dirac field plus the charge due to any external sources in Bp); this is 
easily seen by taking a Gauss integral over the sphere at infinity of the electrostatic 
field (given by the gradient of A°). 

4 No Embedded Eigenvalues 

A famous theorem of H. Weyl asserts the invariance of the essential spectrum of 
the perturbation of an operator if the difference of the resolvents of the perturbed 
and original operators is compact, see §4.3.4 of 0. In standard notation we have, 
for a stationary system with ijj = 



the free operator, Hq, has A“ = 0. Consequently, we have 

Proposition 1 The Dirac Hamiltonian operator H of a stationary and isolated 
Maxwell-Dirac system has the same essential spectrum as the free operator, i.e. 

o'ess{H) = <^ess{Ho) 

= 

= (—oo, —m] U [to, oo). 

Proof A simple adaption of §4.3.4 of Q. □ 

We will now prove that for our stationary isolated systems there are no 
embedded eigenvalues, E (the “energy”). 

Theorem 1 A stationary and isolated Maxwell-Dirac system has no embedded 
eigenvalues E, i.e. —to < E < m. In particular, the following limit exists, 

— = lim —cosy^^ where }? = (F + , 

m kl-oo ^ 

with P = and n“ = 
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Remarks 


• The result needs the Maxwell equations only in order to derive the decay 
result for of lemma |^. The result remains true for the Dirac equation 
alone if we assume the appropriate decay for >1“. 

• The only conditions required are the rather weak regularity and decay con¬ 
ditions of an isolated system. No positivity conditions on the potential are 
required, c/§4.7.2 of 0. 

Proof. The proof is remarkably simple, it is a matter of exploiting the notation 
introduced in definition 0 We begin by re-writing the 2-spinor form of the Dirac 
equations, (0), in this notation. We have, 

+ = 0 ( 11 ) 

^ (~~R~ ^ ~ 

where, = ESq + y4“. 

We combine these two equations into a single (equivalent) equation which 
gives the derivative -I- idaX- To do this multiply the first equation by ts and 
the second by ob and subtract. Using oa^b — ob^a = (ab (see Appendix A), we 
have 


^^R^ + ^^aaX - (i-aob + I'BOa) ( 12 ) 

+2 {iAd\oB - OAd\iB) + V2im {iaIa ~ ^AO^) e~’'^ = 0 . 

We now use the multiplication lemma to place the terms of the equation into an 
appropriate weighted Sobolev space: 

• G m^i^+ 2 e(A'p), since (from definition ||), 

OA G lU|/(Up); 

• (i^os + I'BOa)B^^ = {laob + 1'B0a){E(7^ + A^), breaking this into two 

terms we have, 

{I'AOb + 1'BOa)E<J^ G IU 2 J (Up), {oaOb + I'BOa)^^X ^ ^-’l+2e('^p)’ 

this last inclusion uses lemma Q 

• LaQ^^Ob - OaB^j^Lb G lU^f+2e(-^p)- 
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So that equation implies, 

- ‘^'iii-AOB + iBOA)Ea° \ (13) 

+V2im {iAi^A ~ oaoa) G 

We now contract equation dil) with and so on. Using, the results of Appendix 
A (in particular the last two facts), we have, after splitting the resulting equations 
into real and imaginary parts, 

3 R 

+ V2msinx (ria - la) (14) 

H 

+i{m^ — m°)E{m^fna — rffrua) G W^^j^^^{Ep), 

(Z° — n°), (m° + m°)U, (Z° + n°)E — v^mcosy G W‘^f^^^{Ep). (15) 

Now (using Appendix A), 

1^ + 71^ = --j= (|oo|^ + |oi|^ + |to|^ + ki|^) > 0 and 

3 

(/“ + n“)?7„/3(/^ + n^) = { 1 ° + n°)2 - ^(U + = 2, 

i=i 

here rjafs is the Minkowski metric. So we have. 


■^(Z° + n°) — ■^(^0 + ’^o) — 




l + lU(i'+>^V = l/l + 5V. 


i=i 


From the last inclusion of equation dil) we have 

-^(;° + n°)U — mcosx = 1 + E — mcosx G lU^f^ 4 g(Up). 

Notice that, (1 + ^A^)“5 g and consequently. 


E 

m 


cosy 


G W 


2,2 


Finally, from the Sobolev inequality we have, 


-l+6e 


(Ep). 


E 

m 


cosx 


1 + 


<CU|-^+®^ 


for any e > 0 and some constant C. Hence the limit |a::| ^ oo of the left hand side 
exists and is zero, which completes the proof. □ 
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5 Regularity and Decay 

A stationary Maxwell-Dirac system is an elliptic system of partial differential equa¬ 
tions. So it should be a simple matter to apply the theory of elliptic regularity to 
obtain the best possible regularity results for the Maxwell and Dirac fields. That 
this is indeed the case is demonstrated in the next theorem. 

Theorem 2 A stationary Maxwell-Dirac system for which, Ua, Va G L‘^{Ep) 
and G Lj^^{Ep), is , i.e. U, V and A are in C°^{Ep). 

Proof. First we note that we can always find a gauge transformation which takes 
A into the Lorenz gauge while leaving it in the same Lebesgue space. 

We have a set of Poisson equations for the A’s (equation ®) , and the elliptic 
Klein-Gordon equations (B.3). These equations take the form, 

= 47rej“ 

3 

AUa + 2ie^A>djUA = {m^ - E^)Ua + h^Us, 

j=i 

where h^ is a quadratic function of the A’s; with a similar equation for the P’s. As 

U and V are in L‘^{Ep), is in L^{Ep), so the A’s must be in L^^^^{Ep). Putting 

this information into the U and V equations we conclude (after the use of the 

. 6 

Holder inequality) that the hgllA are in L^^^{Ep) and so the H’s and P’s are in 
L^o(.{Ep) (elliptic regularity). We can now conclude that all the fields are in C°’“ 
and then in (7^’“. Iterating, we finally find that the fields are in C°°{Ep) (this is 
the classical ‘bootstrap’ argument, see chapter 10 of or chapter 9 of |^). □ 
In the case of the free Klein-Gordon equation it is easy to see that, in the 
stationary case, if E"^ — mf < 0 then the Dirac field decays exponentially. This fact 
remains true for our stationary isolated systems. 

Theorem 3 For an isolated and stationary Maxwell-Dirac system, with E^ vrif, 
the Dirac fields U and V (along with all their derivatives) decay exponentially as 
\x\ oo. 

Proof. With E'^ m? and theorem 1 we have \E\ < m. 

We now cast our Klein-Gordan equations in a form suitable for the application 
of the maximum principle. Firstly, A(|[/p) = UAU -\- UAU -\- 2|VC/p, so using 
equation (B.3): 

A(|t/oP + |C/iP) ( 16 ) 


2ieY)^j^i{UoA^djUo - UqA^ djUoUiA^djUi - UiA^djUi) 

+2{m^ - E^){\Uo\^ + |C/iP) + 2(|VC/oP + iVPip) 

-2 [2eEA° + e2A“A„] (|[/oP + |PiP) - ie{dA - dA){\Uo\^ - \Ui\^) 
-2ie [{-d,A + dA^)UoUi - {-d.A + dA^)UoUi] 
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We will use the following comparison function, 


_2 I 1 

w{x) = Cq- ——j— with Co = sup [h(a;)] and k > 0, 

\^\ \^\=p 

where h{x) = (|C/oP + |Cip + |VoP + l^iP); the supremum in the definition of 
Co is well-defined and hnite since we have by the Sobolev inequality, 0 < h{x) < 
on Ep, for p large enough. Note that we have Aw — 2fc^w = 0, on Ep. 
We show that 

A[/i(x) — w(a;)] — 2k‘^[h{x) — w(a;)] > 0. 

The following inequalities will be needed - in each case we have used the Sobolev 
inequality on the A’s (after use of lemma 1). Firstly, 

3 

(UoA^djUo - UoA^djUoUiA^djUi - UiA^djUi) 

i=i 

= R'^[A^djxi\oo\^+ 1011"^) 

j=i 

-I- i{ooA^djOo — ooA^djOo -I- oiA^djoi — oiA^djoi) 

> “j^^doop + |0rp + kop + kiP) 

= -^(ic/ok + ic/ik + ivbk + iiAk), 

where we have used, 

1 = k'^o^l = k*^oo -i- t^oii < —(|oop + |oik -i- kop + kik)- 
Next, in a similar vein, we have 

-2\2eEA^ + e^A°^AJ > 

-^e{^A-^A){\Uo\^ > -e\dA-dA\{\Uo\'^ + \U^\^) 

> -^W + \Uin 
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- 2ie [{-d,A + dA^)UoUi-{-d^A + dA^)UoUi] 

C 4 


> - 2 .-. 


a 


:\UoUi 


> -r^iW+ \u,n 


here the Cj are positive constants. 

With the use of equation ( 0 ) and the inequalities we have, 


2{m^ 


A[(|C/oP + |C/iP)] 

> 

C 2 C 3 


_ _ cu 

\x\ jirp 


W + \u^n 


-^M^) + 2(|Vt/oP + |VC/ip). 

There is, of course, an entirely similar equation for Adding these 

two equations gives, for every k such that 0 < k < V— E"^, 

A [h{x) — w(a:)] — 2k^[h{x) — w(a;)] 

C 


> 


2{m^ -E^ - k^) - 


h{x) + (|v[/oi" + ivc/ii^ + |vK)r + |vwr) 


> 0, for p large enough. 


Applying the maximum principle on Ep we see that the non-negative maximum 
of h(x) — w{x) must occur at infinity or on |a::| = p. However, 

lim {h{x) — r/;(x)] = 0 and 

|x|—^■oo 

[h{x) - < 0. 

We conclude that [h{x) — w{x)] < 0 on Ep, so that 

_/c 13? I 

iC/oP + iC/ip + < Co '' I I ■ 

So the C/’s and f/’s decay exponentially. 

Differentiating the Klein-Gordon equations we can use the same procedure 
to show that the first derivatives decay exponentially. After taking account of 
theorem 2 we can iterate this procedure once we note that the solution for A“ 
can be written as the sum of an harmonic polynomial (of negative degree) and the 
convolution of and the appropriate Green’s function. □ 

This theorem does not deal with decay of solutions in the case \E\ = m. It is 
clear that the solutions in this case need not decay exponentially. The spherically 
symmetric solution of 0 is a case in point, this solution has \E\ = m and R ~ -j^ 
as |a;| ^ 00 . 

The case \E\ = m is in a sense quite unique, as we will see in the next section. 
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6 |£| = m and Asymptotically Static Systems 

In this section we will prove a theorem which neatly ties together the \E\ = m 
case and the concept of an asymptotically static solution. 

In 0 the idea of a static Maxwell-Dirac system was exploited to show that 
if the system was also stationary and isolated then the system was necessarily 
electrically neutral, with \E\ = m. A static Maxwell-Dirac system is one for which 
(in some Lorentz frame) the spatial components of the Dirac current vanish. With 
the Dirac current written as, 

+v^v^) 

= y2i?(Z“-hn“) 

we would require for a static system that = 0, fc = 1, 2, 3. We will now 

generalise this concept to that of an asymptotically static Maxwell-Dirac system. 

Definition 3 A Maxwell-Dirac system will be called asymptotically static, with 
decay rate k and differentiability index s, if l^ -\- € 

Wf'^{Ep) for some p, k > 0 and k = 1,2, 3. 

So an asymptotically static system decays toward a static system, as \x\ oo. The 
unit vector, (in the direction of the current, j°‘) has only a time-like 

component in the limit as |a:| —> oo. 

To take full advantage of this definition we will need to recast it in terms of 
the individual variables oa, and i-a, this is done in the following lemma. 

Lemma 2 An isolated and stationary Maxwell-Dirac system is asymptotically 
static if and only if 

t°-do, -di&W!:l+^{Ep). 

where t = mm[3, s]. 

Proof. First assume the system is asymptotically static with decay rate k and 
differentiability s 

From the proof of theorem 1 and definition ^ we have, 

(|oo|' + |oi|' + |t°|' + kY)'-4 = 2{f + n°f-4 

= 2^il>^ + n'^rGWlt{Ep). 

k^l 

Consequently, as [(|oop -I- |oip + -\- |t^P) -I- 2]“^ G W 2 ^{Ep), we have 

loop + |oip + |.°p + - 2 G 
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Writing out the explicitly we have, 

— (ooOi + OiOq) + € Wf_’^(Ep) 

— (—ooOi + oiOg) + (— € Wl’^(Ep) 

-|oop + |oip + |,op-|,i|2 e Wl’^(Ep). 

So we conclude that, 

—ooOi + G Wl’^(Ep) and 
|oop + |.i|2-l, \o^\^ + \,^\^-l^whl{Ep). 

Writing these equations as a single matrix equation. 


Ol 

\ 

/ Og-tO 

Oi — b^ 

Oo 


1 

o 

rt 7^ 

Oo — ^ 


oiOp - |oip + |t°p 7 1 \ 

|ooP + k^P-l OQd^-b^lP j 


is in W^J^{Ep). The first matrix on the left has determinant 1 and inverse, 



which is in W^'-^{Ep). The result now follows from the multiplication lemma after 
applying this inverse matrix to the previous equation. 

Next, assuming — dg, — dj G we easily find that + n’^ € 

W^‘^^ 2 ei^p)t k = 1)2,3, and any e > 0. So the system is asymptotically static. □ 
Now to our theorem connecting the two apparently unrelated notions, the 
condition |i?| = m and the idea of an asymptotically static system. 

Theorem 4 A stationary and isolated Maxwell-Dirac system is asymptotically 
static if \E\ = m. 

Proof. Assume the system is stationary and isolated with \E\ = m; write E = em, 
with e = ±1. 

The proof is very simple, it simply involves manipulating expressions obtained 
in the proof of theorem 1. From that proof we have, 

-ecosx G VF!f+ 4 ,(F;p). 

Note that we must have ecosy > 0, for p large enough, so 

+ ecosx G Wl'^{Ep). 

Multiplying the last two expressions and using the multiplication lemma, we have 

sin^X ^ lb-i+8,(-Ep). 
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These equations together with the fact that both sin^ x and l^+n^ are in W|j^(i5p) 
enable us to conclude that 

sinx, eW^^jEp). 

The system is, according to definition 3, asymptotically static. □ 


7 Exponential Decay, the \E\ = m Case 


In this section we will prove that the Dirac held decays exponentially in the |ii^| = 
m case as well - at least when the total charge go 7^ 0- In fact we obtain tight 
bounds on the decay of the Dirac held in this case. 

Theorem 5 The Dirac field for a stationary, isolated and asymptotically static 
Maxwell- Dirac system, with \E\ = m, k > 1, s = 3 (definition 3) and qg ^ 0 
decays exponentially as \x\ oo. In fact, there exists two positive constants Ci 
and C 2 such that 


-4V2mAy^ 

Cl - 3 - <R< C2 — 


|a;|2 


|a;|2 


where A > 0, = —ee— is necessarily positive and — = e = ±1. 

■ rr? a r 777 , 


Before embarking on a proof of this theorem we will need a couple of preparatory 
lemmas. In the course of proving the second of these two lemmas we will also show 
incidentally that the electric dipole moment must vanish. We are assuming that 
\E\ = m and write -^ = e. As cosy ^ e, with |a;| ^ 00 we can take y = utt + ^ 
where (—1)"' = e and C ^ 0, as |a;| ^ 00 . 

We require a more careful analysis of equation (12). The imaginary part of 
this equation is, 


daX +2{E + eA°){n°la - l^Ua) + A^{l^na - n^la) 

+y2mcosy (uq, - la) - i {o^d^^OA - OAd^^dx)na 
+{EdAAE - Ed^MEla + {OaQa^E - L^dA^ojfjma 


(17) 


+ iEd^lOA- OAd^E) 


rUo 


= 0 . 


Assuming the system is asymptotically static and making use of lemma 2 simple 
calculation reveals, 


-Cd.iE + OAd^^o 


- OaOa^E 


'A ^ ^-K-l+2eEp) 


G WE,_,^2eiEp)- 


From lemma 2, G , 2 j(£'p). 
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For a = A: = 1, 2, 3 we make use of G Wl’^{Ep) and our previous 

equation to get 

2{n°lk - fuk) - {1° + n°){lk - Uk) = (n° - f){lk + Uk) G W^t.+eAEp). 
Using lemma 1 and the asymptotic staticity of the system 

3 

2eY,AHPnk - nHk) e W^jl_,^^,{Ep). 
i=i 

Putting this all together in the equation for dkX we have 


9fcC + 


{em + + n°) — emcos^ 


{Ik klk) G VF_’i_„_|_gj(£'p)- 


We need to refine our estimate for (C — n^)l\/2, noting 


fl° + n° 


1 ^ 


■n^ f e W^_t{Ep), 


and [—^{l^ + u°) + 1] ^ G W^fiEp) we have 

v2 




Now use this estimate, together with lemma 1 (to separate the monopole term) 
and the fact that {Ik — rik) + 21^ G w!^"^{Ep) (remember, Ik = —l^), to get 


dkC — 2'j2em 


1 — cos + 


gego 

m\x\ 



m 


l'^ G W 


2,2 

-i^+9e 


{Ep), 


where v = min[2K, 1 + k]. For convenience we will write this equation in a 3-vector 
notation, using V to denote the gradient and I = {P,P,P); note that (because 
of asymptotic staticity) we have for the norm of I, |Zp =1.1 = ^ plus a term in 
(see proof of lemma 2). We have. 


VC - 2V2f 


em 


1 — cos C + 


eeqo ee 

—n ^ 

m\x\ m 


I G Wl’^,,,{Ep). 


(18) 


We will also require the asymptotically static version of equation (^, the 
real part of equation (ED- Firstly, from lemma 2 we have that 


V2m^ = V2a^^j^o^L 

= — OqI^ + Oit® 

G Wl’^^^,{Ep). 
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We have, after using, {Ik — nu) + “21^ G W^'^{Ep), 

ivi? + 2v^emsinCZ S (19) 

Now, to the first of our two lemmas. In fact, this lemma actually gives us the 
exponential decay result which we use to improve lemma 1 so that we may obtain 
the much tighter estimate necessary for theorem ||. 

Lemma 3 The Dirac field (and at least its first and second derivatives) of a sta¬ 
tionary, isolated, asymptotically static (with k > 1, s = 3 and \E\ = m) Maxwell- 
Dirac field decays exponentially as |a;| ^ oo, provided qo 0. In particular. 


g-2fc-\/H 

R < C -;—;-, for any k such that 0 < k < 2v2mA, 

\x\ 

here is necessarily positive. 

Proof We first note from theorem 4 that sin ( G so, for p large enough, 

C G W"!|+4e(^p) and VC G W]:l^^^{Ep). 

Now, taking the divergence of equation and using (H) again to remove 
the |Vi?p term. 


— Ai? — Avnf sin^ C G W^'s {Ep) for any e > 0. 
From the fact that VC G W]_''l_^^^{Ep) we have from 


1 — cos C + 


eeqo ee q _ . 2 C , eego , ee g 


0=2 sin' 
m\x\ m 2 


777 X m 


+ G 


Our first observation is that eeqo < Oj since from lemma 1 a° S with rj > 1. 
We write and take A > 0. We can also use the last inclusion to 

m ’ 

estimate the term sin^ C- We have, 

sin^C - ^ G \ .3 1, (Ep). 

The second order elliptic equation for R can now be written as 


AR -8m^ — 


f + o') i? = 0, where a G . .3 , (Ep). 


We will now use the maximum principle utilising a comparison function 

k 


g-2fc-\/N (^2 

v(x) = C -j—j-, for which Av — 


X 


X 


2|a:|i 


u = 0. 


(20) 
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Now, 


A[R — v(x)] — 8m^ 


+ a ) [R- t>(a:)] = 


(8m^X^ — fc^)— + a 
\x\ 


}[x), 


where a € r)]+e^^P^' Consequently, for p large enough and for every k 

such that 0 < A: < 2^/2m\ we have. 


A[R — u(a;)] — 8m^ 



[i? — u(a;)] > 0. 


Choosing C such that [R — v{x)]\x\^p < 0 we have by the maximum principle that 
R — v(x) < 0 on Ep. Completing the proof of our lemma. □ 


Lemma 4 For a stationary, isolated asymptotically static Maxwell-Dirac system 
with K > 1, s = 3 and \E\ = m we have the following estimates when qo ^ 0. 


£iV 2A £i £i( 16AW + 9) 3^2 .p X 

4m|x| 96v^Am2|x|i ^ 

= (f — u), with u € bCf f 


rd -^ e and f.u G 

v2 

where f is the radial unit vector and (£i)^ = 1. 


Proof In all that follows we are assuming that p is large enough that the necessary 
expansions - eg sin ^ G , (^p) when sin^ G {Ep) - can be made 

on Ep. 

We begin with the estimate. 


sin^C- ^ e ,3 , 

3, -min[f ,r;]+e 


from the proof of lemma 0. Write, 


_ V^£^ ^ where £1 = ±1, 


and substitute into equation 


r X 

VCi- ^J—^f — 2V2em 


V2| 


x \2 


2V2S1X 

'JW\ 


Cl + Cl" 


I G W%,iEp). 
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We have kept only terms “less than order l/jxp” on the left of the equation. The 
from is of order l/lxp, since we can now improve the result of lemma 1 using 
lemma || - from equation we have = 47re-\/2(Z° + n°)i? so must be 

the sum of an harmonic polynomial (of negative degree) and a term which decays 
exponentially. 

Starting with (^i G (since ( is) we have that VCi G W ^’3 , in the 

2 2 

first instance. But then our equation (above) implies that Ci G and that, 

jijf + 4:V2m£(i I G 
\x\ 

Now write Ci = ~ 4 m\x\ ‘^2 {^2 = if) and repeat the process to find that (2 € 

W^’l . As a consequence we have the following estimates for ( and I, 

2 


c- 


££2 


V2X£i 
\/R 4 m|a;| 

£2 


e W°i(Ep) 




I = —= {f — u), with u G W ’t. {Ep). 
3/2 " 2 +" 


Next we use the fact that, LI — i G as the system is asymptotically static. 

We have, 

3202 32 

It is a simple matter to show that for u G W 2 ^ n we have u G W_\ 

Begin with a function / in n then integrate V. (o'^'^“^/V/) over 

Ep to show that V/ G for 0 < d < iterate this process to eventually 

find V/ G W^’l which together with / G shows that / G W]^1 . 

Now repeat the process with dkf G 1 ^ fi , , and so on to eventually get 
/ G W^'i (Ep). With u G , we can use the multiplication lemma to get 

|mP = u.u G .Next we use the fact that, /.Z — 4 G as the system is 

asymptotically static. We have, 

l.l — - = -{—2f.u + |u|^) G 


So that f.u G as K > 1. 

We also note that an argument similar to that used above shows that as 


C2 G n 


nw 


0,2 


-4+e 


3 2 

we must have C 2 G WJ 3 . We can now substitute, 

9 


C2 = + Cs, where ao = £1 

a; 2 


(IGA^to^ + 9 — 96m^d.f) 
96V2Am2 
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here d is a constant vector arising from the expansion of = -[^ + ]^+ 0 (y^). 

We find that (3 G We note from equation ( p^ that 0.V( and 0.VC are both 

in W^ 2 _f_^ - here 0 and are the angular unit vectors orthogonal to f. Consequently 
we must have d = 0 . 

So, as d gives rise to the electric dipole moment, we see that the electric dipole 
moment must vanish. 

We still have to show that £2 = ££i to obtain the precise statements of the 
lemma. This is easily done by taking the estimates for sin ( and I and substituting 
them into equation ( 0 ) - to obtain exponential decay (rather than growth!) we 
require ££162 = 1 . 

□ 

Proof of Theorem 5 Armed with lemma ^ the theorem is remarkably simple to 
prove. We start with equation (||), which can be written as 

—hVi? + v^ETOsinC + V.^ = 0. 


Which gives, 


r.Vi? M.Vi? „ .. 2 ^ 

-h 2£imsm4 H-V.n = 0. 

R R r 


Using equation 0 ) we have (as u.l G 


3.2 N 
l+e) 


u.VR 

R 


G W 



Noting that V.u G , we have, using lemma 0, 

2 


f.VR 

R 


„ /- mX 3 


G lU^’a%,(Ap). 


Consequently, we have 


R = K- 


„-4V2mAy^+/3 


where [3 G and K may depend on From equation ( p^ we see that 

IniF G so that AT is a bounded function. The Sobolev inequality implies 

\(3\ < ( 7 /|a;|“ 2 +'^. The result now follows by bounding, Ke^. □ 


8 Discussion 

It is worth emphasising here that our results are based purely on rather weak 
asymptotic regularity and decay assumptions. Nothing is assumed about the be¬ 
haviour of the fields in the interior region Bp. 
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Another important point to note is that all the results require the Maxwell 
equations only to obtain the decay conditions of lemma 1 and the improved decay 
required for theorem |^. If this decay is given, a priori, then the results apply to 
the “Dirac equation in an external field” as it is usually presented. 

One question which needs to be addressed is the possible extension of the 
electric neutrality theorem of to the asymptotically static case. But we leave 
this to a future paper. 


Appendix A: 2-Spinors and Spinor Dyads 

We collect here a number of facts relating to 2-spinor dyads and their associated 
null vectors. We give only a brief statement of the facts, for details the reader 
should consult the book of Penrose and Rindler, |l^ . 

• 2-spinor indices are raised and lowered with and cab (summation on 

repeated indices), and ^a = (-ba^^ for any 2-spinor . 

T A B BA AB ^ 

• oa^b — OBiA = ^AB and o L — o i = e , where 

(cab) = (e^^) = ( J ) • 

• OA^^ = —o"^LA = 1 and = 1. 

• The van der Waerden symbols connect Minkowski vectors to 2-spinors 

and vice versa. The \/2{aA^) (with k = 1,2,3) are simply the Pauli ma¬ 
trices and \/2{aA^) is the identity matrix. We have, = pap the 

Minkowski metric, and aA'^°‘bb ~ 

• Because of these relations the null vectors P = cr“ n“ = cr“ 

m“ = and m“ = form a null tetrad; with la and ria real, 

and TTia complex. We have, l°‘la = 0, rf^ria = 0, m°‘ma = 0, PtUq = 0, 
l°^fha = 0, n°‘ma = 0, n°^fha = 0, l°‘na = 1 and m°^fha = —1. 

• For any vector Xa we have, 

Xa = {n^Xp)la + {l^Xp)na - {fh^Xp)ma - {m^Xp)fha. 

• = rdaf, = n°‘daf, and so on. 
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Appendix B: Explicit Forms of the Dirac Equations 

In this appendix we collect together explicit forms of the Dirac and Klein-Gordon 
equations for stationary Maxwell-Dirac systems. 

d 

In this section we use the notation 9^ = — etc ,d = dx + idy, and for the 

oz 

electromagnetic potential A = + iA^. 

The Dirac bi-spinor is, 




Ua \ 

)- 


The Dirac equations are, 
E 


im{y° - —UA - dUi - d^Uo - ie -b AAUq + AUi] = 0 
m 

E 
im 


(BA) 


(yi _ ^Ui) - dUo + d,Ui - ie [AUq + - A^)Ui] = 

m. 


0 


m 

im{Ui - —V^) + dV° - + ie 

m 


-{A° + A^)V'^ + AV^ 


= 0 
= 0 


im{Uo -y°) -b dV^ + + ie AV^ - (A° - AAV° 

m L 

The Klein-Gordon equations are easily derived via differentiation of equations 
(B.l), we give the results for U only, 

AC7o + AdjUo + {{E^ - w?) -b 2e£;^° -b e^A°‘Ao, (B.2) 

+ ie [d,{A° + A^) + BA] }Uo + ie [d,A + d{A° - A^)] Ui = 0 


AUi + AdjUi + {{E'^- rrA) -b 2eEA^ + 

-b le [-d,{A° - A^) + BA] }Ui+ie [-B^A -b 5(^0 -b C/q = 0 

Equation (^) gives the electromagnetic potential in terms of the C/’s and K’s 
which may in turn be written terms of R, x the o’s and t’s. We give only the 
result for A°, 


m [g-zX(|oo|2 + |oi|2 + |,0|2 + |,l|2)_ 


2E1 


(BA) 


+ iBxVoi + B{Aoi) + (^ + iBx)l^oo + B{Aoo) 

+(^ + + Bz{l°oo) - -b iBzx)‘'^oi - Bz^b^oi)] . 
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